5 = 1/2 chains and spin-Peierls transition in TiOCl 
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We study TiOCl as an example of an S — 1/2 layered Mott insulator. From our analysis of new 
susceptibility data, combined with LDA and LDA+U band structure calculations, we conclude that 
orbital ordering produces quasi-one-dimensional spin chains and that TiOCl is a new example of 
Heisenberg-chains which undergo a spin-Peierls transition. The energy scale is an order of magnitude 
larger than that of previously known examples. The effects of non-magnetic Sc impurities are 
explained using a model of broken finite chains. 
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The discovery of high-T c superconductors has gener- 
ated a great deal of interest in low-dimensional spin- 
1/2 materials. Even in the absence of charge degrees 
of freedom, many questions remain unanswered until to- 
day. One of key importance is the question whether a 
two-dimensional spin system can support unbroken sym- 
metry at zero temperature due to strong quantum fluc- 
tuations for S = 1/2 systems. Such a spin liquid state 
was suggested by Andersoro and termed the resonating 
valence bond (RVB) state . This has motivated us to 
search for other examples of S = 1/2 layered materials, 
notably at the beginning of the transition metal series 
where Ti 3+ is in the d 1 configuration (as opposed to d 9 
in the cuprates). From this point of view the layered 
compound titanium oxihalides TiOX (X=C1, Br) appear 
most promising. Indeed, Beynon and Wilsona reported 
in 1993 that the uniform magnetic susceptibility of these 
materials shows a number of unusual properties. Accord- 
ing to them, the susceptibility is almost temperature in- 
dependent and could not be fitted to the Curie- Weiss 
law, nor to any one-dimensional model. Moreover, if 
non-magnetic Sc (d°) impurities were introduced, a large 
Curie tail corresponding to one spin- 1/2 moment per Sc 
appears. After subtracting this Curie tail, the suscep- 
tibility was again found to be temperature independent 
but substantially reduced compared to the pure mate- 
rial. This effect was particularly striking in TiOCl and 
led Beynon and Wilson to propose that these materials 
may be examples of RVB-like states. 

We have redone the measurement of the susceptibility 
of TiOCl for both the pure material and in the presence 
of Sc, and carried out LDA and LDA+U band structure 
calculations. Our data are not in agreement with that 
of ref.|[ especially for the pure sample. Here, we find 
that above a temperature of 130if , below which the sus- 
ceptibility drops abruptly, our data fits well to a nearest 
neighbor Heisenberg model with an exchange constant of 
J = 6WK. The data for 6% and 10% Sc-doping is in 
very good agreement with a model based on broken fi- 
nite linear Heisenberg-chains. 

In the following we describe the structure of TiOCl and 



discuss the possibility of Id spin chains formed by t2 g or- 
bitals. The reasoning here will be justified by band struc- 
ture calculations, which feature isolated one-dimensional 
(i-bands. We then present our experimental data for the 
pure and Sc-doped samples and fits to theoretical mod- 
els. We will argue that all evidence is consistent with a 
picture based on quasi-one-dimensional spin-1/2 chains, 
where the pure material undergoes a spin-Peierls transi- 
tion at low temperatures. 

The structure of TiOCl is of the FeOCl-type which 
consists of bilayers as displayed in Fig.0a). The symme- 
try is orthorhombic with Z = 2. The layers repeat in 
the c-direction of the crystal with c = 8.03A. CI layers 
mediate a weak van der Waals interaction between suc- 
cessive bilayers. Within each bilayer, Ti and O form two 
layers of buckled chains where Ti is always on the outer 
side with respect to the bilayer. The O-Ti-0 bond angle 
is 153°. Each Ti ion is surrounded by a distorted oc- 
tahedron of O and CI ions (Fig.[l]a))). These octahedra 
are formed by two O ions belonging to the same Ti-0 
chain, two O ions belonging to neighboring chains, and 
two CI which lie on the outside of the bilayer. They are 
corner-sharing in the a-direction along the Ti-0 chains 
and edge-sharing in the b-direction. Figs.|l]b) and |l|c) 
show a top view of the lattice showing only the Ti ions 
for clarity. The Ti sublattice consists of two rectangular 
layers with lattice parameters a = 3. 79 A and b = 3. 38 A 
in each layer. The top layer is shifted laterally and is 
displaced vertically from the bottom layer by 1.96 A. The 
shortest Ti-Ti bond length turns out to be the distance 
between Ti in different layers, and at 3.21A is just slightly 
shorter than the Ti-Ti distance along b. For orientation, 
the Ti ions in each layer are bridged by oxygen ions (not 
shown in figs.[l]b)+c)) along the a-direction, forming the 
buckled Ti-0 chains described earlier. Note that these 
buckled chains (Fig.^Ja)) are introduced for the purpose 
of describing the crystal structure only, and should not 
be confused with the spin chains, which we will next pro- 
pose, that are based on the electronic structure. In par- 
ticular, we believe that the important exchange path is 
direct t2 g orbital overlap, rather than superexchange 
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FIG. 1: a) Structure of TiOCl. A single bilayer is shown. Note that each Ti is surrounded by 4 O and 2 CI, forming a distorted 
octahedron. These octahedra are corner-sharing along a and edge-sharing along b. b+c) Top view on the lattice showing only 
Ti ions which form two rectangular layers. The smaller symbols are for the bottom layer, b) d xy orbitals forming linear chains, 
c) d xz orbitals forming zig-zag chains. Note that the plane of the orbital is tilted by 45° out of the paper and two of its lobes 



point directly at the Ti ions on the adjacent layer, 
ordering will be needed to break the symmetry. 



Another zig-zag chain is formed by the degenerate d yz orbital and orbital 
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FIG. 2: LDA-DOS projected onto Ti, O and CI orbitals 

via oxygen. At each Ti site, the approximately octahe- 
dral co-ordination dictates a set of axes for the conven- 
tional e g and ti g orbitals. We chose z = a, and x and 
y axes which are rotated by 45° relative to the b and c 
axes. With this choice of axes, the d x i_ y i and d z i or- 
bitals point towards the oxygen and chlorine neighbors 
as required for e g orbitals. The important orbitals are 
the t2 g 1 s (d xy , d xz and d yz ). The d xy orbitals form lin- 
ear chains running along the b-axis, linking Ti ions in 
the same layer. This is shown in Fig.[l|b). On the other 
hand, Fig|yc) shows that the d xz orbitals are tilted by 
45° out of the plane and two of the lobes point directly 
at the Ti ions in the other layer. This is also the direction 
of the shortest Ti-Ti distance and we propose that these 
orbitals may form a zig-zag chain along the a-direction. 
Note that such a state is degenerate with a similar state 
derived from the d yz - orbitals, where the latter are con- 
necting a different set of pairs of neighboring chains. A 
state with infinite zig-zag chains is thus one of a broken 
two- fold symmetry between the d xz - and d yz -orbitals. 

The arguments given thus far assume that the states 
at the Fermi level are ^-like electrons. In order to con- 
firm this, and to investigate orbited ordering, we have 
performed both LDA and LDA+UB calculations using 
the full-potential LMTOu method. j-All calculations were 
performed using the LmtART codeD, and the experimen- 
tally determined unit celfl, and ferromagnetic spin polar- 
ization. Both LDA and LDA+U calculations yielded a 




FIG. 3: a) LDA band structure calculation for TiOCl. 
The Brillouin zone is orthorhombic with X — 2tvx /a,Y — 
27vy/b,Z — 2tyz/c and S is the zone corner. The unit cell con- 
tains 2 Ti atoms. 6 bands around the Fermi level are identified 
as t2g bands well separated from the e g bands, b) LDA+U 
band structure calculation with split off d xy bands. 

magnetic moment of per formula unit. In the follow- 
ing discussion, the band structures and DOS correspond 
to the majority spin, as the minority spin for the d-states 
are completely unoccupied. We begin by presenting the 
site-decomposed density of states Fig.^. As shown, the 
oxygen and chlorine p-levels form well separated bands 
from the Ti d-levels with only small hybridization be- 
tween the two. Also, the octahedral crystal field has 
clearly split the d-states into t<i g and e g contributions, 
with the Fermi level lying within the ti g peak as expected 
for a d 1 configuration. These features can also be seen in 
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FIG. 4: Susceptibility of crushed crystals of pure TiOCl. The 
inset shows single crystal data with and without Curie sub- 
traction. 

the band structure (Fig.^ja)). Therefore, the assumption 
of t2^-like electrons at the Fermi level is valid. However, 
LDA predicts this material to be metallic, and project- 
ing the DOS onto the ti g orbitals indicates that all three 
orbitals are partially occupied (not shown). In order to 
go beyond the LDA and explicitly include the effect of 
strong on-site interactions which may induce orbital or- 
dering, we have carried out LDA+U calculations. We 
performed calculations with U = 3.3eV and a ferromag- 
netic on-site exchange of leV (Figj^b)). As shown, two 
nearly degenerate, one-dimensional bands split off from 
the rest of the t<i g bands creating an insulator (note that 
there are 2 atoms per unit cell). These two bands are 
derived from the d xy orbitals corresponding to the linear 
chains in Fig.|l|b). Varying U only had a significant effect 
on the splitting between the occupied and unoccupied 
bands, not on the shape or width of the occupied bands. 
If the band width of 0.9 eV is identified with 4t, where t is 
the nearest neighbor hopping in a one-dimensional tight 
binding model, we can estimate the exchange constant in 
an effective Heisenberg model via J = At 2 jJJ = 714K. 

Other band-structure calculations for TiOCl found in 
the literatureLl that are based on a phenomenological 
tight binding model are not in agreement with our LDA 
calculations. In ref.[?], CI seems to play an important role 
at the Fermi surface and strong correlation effects leading 
to orbital ordering have not been included. 

The band structure thus indeed suggests the existence 
of one-dimensional spin-1/2 chains in TiOCl, where the 
spins are localized in Ti ^-orbitals. We have repeated 
the measurement of the uniform susceptibility of pure Ti- 
OCl (Fig.[|) and in the presence of Sc impurities (FigJ|), 
over a temperature range from a few K to 800K. 

Single crystals TiOCl and Tii-^Sc^OCl were prepared 
by the chemical vapor transport method using excess 
TiCl3 as transport agentB. TiOCl precursor is a mixture 
of Ti02/TiCl3 (1:6 by weight) and Sc doped precursor is 
a mixture of TisOa/ScsOa/TiCla with 6% and 15% of Sc 
in the evacuated sealed quartz tubes. TiOCl crystal is 



FIG. 5: Susceptibility of Sc doped TiOCl (batch A) 

transported by TiCl3 vapor to the cold zone with a gra- 
dient maintained at 650C/550C within 20cm. It takes 
three to five days to complete the vapor transportation. 
The Sc doping level of the single crystal may not follow 
the calculated initial Sc/Ti ratio due to the uncontrolled 
transport mechanism with excess TiCl3. However, we 
find that the thoroughly mixed and reacted powder sam- 
ples have Curie constants corresponding to the amount 
of free Ti 3+ spins that are very close to half of the calcu- 
lated initial Sc percentage. 

The average Curie constants for single crystals are 5.8± 
1.2% Ti 3+ for 15% Sc initial mixture (batch A) and 3.1 ± 
0.3% for 6% initial Sc mixture (batch B). It is very likely 
that the Curie constant does reflect only half of the Sc 
doped in the sample. 

Our data for pure TiOCl shown in Fig.|| were obtained 
from crushed crystals. They are about a factor 3 smaller 
than those proposed in ref.||. After subtraction for a 
small Curie tail and correction for trace amounts of fer- 
romagnetism, we find a sharp drop in the susceptibil- 
ity at low temperatures. In a picture based on Id spin 
chains, such a drop must be interpreted as a spin-Peierls 
transition. The fiat part of the curve below 50K then 
determines the zero-level of susceptibility of the Ti d- 
electrons. The single crystal data shown in the insert 
of Fig. [I] features a sudden drop of the susceptibility to 
zero at 67 K. In addition, a noticealbe inflection point 
appears at 95 K. The high temperature part of the data 
(above 130K) can be fitted to the Bonner-Fisher-curvdj, 
using the nearest neighbor exchange J as the only free 
parameter of the Heisenberg-Hamiltonian: 

H = jJ2Si-S i+1 (1) 

i 

where the Si are spin-1/2 operators. The solid curve 
in Fig.|| shows that the data are indeed very well de- 
scribed by this model. Note that both the absolute mag- 
nitude and the peak position are well accounted for by 
a single parameter J, justifying the choice of zero made 
above. The fit determines J = 660if, which agrees well 
with the crude estimate given above. The result of Sc- 
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doping can now be explained. In our model, we assume 
that Sc replaces Ti on a fraction x of all Ti-sites. Sc 
is in a d° configuration, therefore Sc-doping can be re- 
garded as the analogue of Zn-doping in the cuprates. The 
Heisenberg-chains of our model are thus broken into fi- 
nite open-ended chains of an average length of 1/x. More 
precisely, let P x (n) = nx 2 (l — x) n denote the probability 
for a site to belong to a chain of length n. Denoting the 
susceptibility per site of a finite open-ended Heisenberg 
chain of length n by x( n iT), we expect the measured 
susceptibility in the presence of Sc to be of the form: 

oo 

X sc(x,T) = Y, p *(n)x(n,T) (2) 

n=l 

We have determined x( n ^) by exact diagonalization 
for n = 1..16 and extrapolation. As shown in Fig. ||, 
this model with x = .1 agrees well with the experiment 
on batch A at all temperatures below 600i^, above 
which the doped sample begins to decompose. Similar 
agreement was achieved for batch B with x=.06. The 
existence of a Curie-like tail at low temperatures is easily 
understood in our model and is due to the presence of 
chains of odd length, which behave as free spin-l/2's at 
temperatures below the gap. Since the gap decreases as 
1/n as the chain length n increases, chains of different 
lengths will contribute to the Curie-like tail at different 
temperatures. The fitting to a Curie-tail is therefore 
imperfect, as pointed out in ref.||, except at very low 
temperatures, where all odd chains contribute. It is 
clear that in the limit of small x half of all chains in 
the system will be odd. Hence in our model, Sc doping 
roughly introduces free spin-l/2's into the system on 
a 2:1 basis at very low temperatures, as conjectured 
before. This result is in contradiction with refj2| where 
a 1:1 correspondence is reported. Furthermore, in ref.|2| 
a large drop of the flat high temperature part of the 
susceptibility was found upon introduction of Sc into 
TiOCl. This again does not agree with our measure- 
ment, where the high temperature parts of the pure 
and doped sample have very similar values over a wide 
temperature regime. We remark however that the effect 
of Sc doing as measured by us resembles that reported 
in ref.|| for the case of TiOBr, as one would expect 
for two such similar chemical compounds. Note that 
our model does not assume any spin-Peierls transition 
for the Sc-doped cases. However, we believe that the 
susceptibility data alone does not allow us to draw any 
firm conclusion as to whether the spin-Peierls transition 
is destroyed or not. 

We conclude that both our numerical and experi- 
mental results consistently indicate the existence of 
one-dimensional spin-1/2 chains in TiOCl. Above 130if , 



the susceptibility of TiOCl is well described by a nearest- 
neighbor-exchange Heisenberg model with J = 660K. 
We interpret the abrupt drop of the susceptibility at 
67K as a transition into a spin-Peierls state. We have 
identified the two possible kinds of spin chains that may 
be present. These are linear chains along the b-axis and 
zig-zag chains along the a-axis of the crystal. The ti g 
orbitals that form these chains have been shown to be 
energetically isolated in the band structure calculation. 
We note that in the case of the linear chains it is clear 
from fig.|l|b) that the neighboring chains form a staggered 
pattern and the exchange coupling between neighboring 
chains is frustrated. Thus interchain coupling effects are 
likely to be weak. Furthermore, below the spin-Peierls 
transition, the lattice distortion associated with the 
dimerization is also frustrated from one chain to the 
next. This is not the case for the zig-zag chains. Once 
the degeneracy of the zig-zag direction is broken, the 
interchain exchange is unfrustrated and zig-zag chains 
can easily lock in place. Thus interchain effects are 
likely to be different in the two scenarios. A first order 
phase transition into a Spin-Peierls state as suggested by 
our single crystal data is expected by Ginzburg-Landau 
arguments only in the frustrated scenario of the linear 
chains. This is also the scenario favored by the LDA+U 
calculation. Further experiments such as x-ray diffration 
are needed to confirm this picture. 

A number of other spin-Peierls systems have beep, 
studied, notably the organic system TTF-CjjBDTlI 
(where the spin resides on TTF) and CuGeOslij. The 
exchange constants and transition temperatures are 
J = 77K, T c = 12K for TTF-CuBDT and J = SSK, 
T c = 1AK for CuGe03 and the interchain coupling is 
unfrustrated in both cases. The energy scale for TiOCl 
is close to an order of magnitude higher. The relatively 
large J and the weak Jahn- Teller effect for t^g orbitals 
both tend to minimize carrier localization effects making 
TiOCl an attractive candidate for doping. For large dop- 
ing x where xt exceeds the spin-Peierls energy scale, the 
doped holes will destroy the spin-Peierls state, resulting 
in a Luttinger liquid. For smaller hole concentration the 
spin-Peierls local order survives and the spins remain 
in singlet pairs. Coherent motion of doped carriers 
may lead to unconventional superconductivity, as in the 
original RVB scenario. These possibilities are certainly 
worthy of further investigations. 
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